Superfluid-insulator transition and BCS-BEC crossover in "dirty" ultracold Fermi gas 
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Superfluid-insulator transition in an ultracold Fermi gas in the external disorder potential of the 
amplitude Vo is studied as a function of the concentration of the gas n and magnetic field B in 
the presence of the Feshbach resonance. We flnd the zero temperature phase diagrams in the plane 
{B, n) at a given Vo and in the plane {Vo, n) at a given B. Our results for BEG side of the diagram 
are also valid for the superfluid-insulator transition in a Bose gas. 



Using the Feshbach resonance in the magnetic field 
B one can study lots of interesting physics in ultracold 
Fermi gases (see the recent review article^ and refer- 
ences therein). In the vicinity of the Feshbach resonance 
B = Bq the scattering length of two fermions typically 
changes as 



a = ao 



AB 

Bn-B' 



(1) 



where we omitted the non-resonant term. As a result 
by the decreasing magnetic field the Fermi gas can be 
transformed from the phase of weakly attracting fermions 
(at B > Bo) to the phase of repelling each other com- 
pact composite bosons, dimers made of two fermions 
with opposite spins {aX B < Bq). At B = Bq the gas 
goes through unitarityi. In a clean Fermi gas all men- 
tioned above phases are superfluid. Far enough from the 
resonance aX B > Bq superfluidity is described by the 
Bardeen-Cooper-Schricffer (ECS) theory, while on the 
other side, at B < Bo, the theory of Bose-Einstcin con- 
densation (BEC) of composite bosons works. Thus, re- 
duction of magnetic field B leads the gas through BCS- 
BEC crossover. 

The aim of this paper is to consider the zero temper- 
ature BCS-BEC crossover in a "dirty" Fermi gas, i.e. 
in the gas situated in a three-dimensional (3D) random 
potential. Such a random potential can be created, for 
example, by superposing a 3D speckle on the ultracold 
gas sitting in a trap. Obviously, a strong enough dis- 
order can localize the Fermi gas on BCS side and the 
BEC condensate on BEC side, destroying superfluidity 
in both cases. For brevity, we call the localized phase 
"insulator" and the localization transition "superfluid- 
insulator" (SI) transition. In this paper we are talking 
about SI transition in uniform in average infinite gas but 
some of our results can be applicable to experiments with 
wide enough traps. 

Expansion of BEC condensate of ultracold Bose gases 
in the disorder potential of one-dimensional speckles has 
been recently studied experimentally^!^. It was found 
that the disorder stops expansion at some distance. In 
this case, however, a big role may be played by rare very 
high hills of the random potential. Apparently several 
laboratories are planning similar studies of SI transition 
in a potential created by 3D speckles. One can expect 
that in this case the rare high hills are less important and 
theory of SI transition in infinite system is more relevant. 



In a Fermi gas in a fixed external random potential 
the SI transition can be driven by the decreasing concen- 
tration of fermions n at a given magnetic field B, or by 
the decreasing _B at a given n. Therefore, one can think 
about the SI phase diagram of a Fermi gas in the plane 
{B,n). In this paper, we find the zero temperature SI 
border line n{B) on such a phase diagram (see Fig. [T]). 
Because B and n can be independently controlled this 
diagram can be verified experimentally. We character- 
ize disorder by the amplitude of the random potential 
energy Vq (mean square deviation of random potential 
V{r) from average value) and the characteristic size of 
potential wells and hills R. In the first part of this paper 
we assume that both Vq and R are so large that if m is 
the mass of the fermion 



mR"^ 



(2) 



This means that a typical potential well has many levels. 
In this sense we are talking about classical random poten- 
tial. From the beginning we assume that n{B) is so large 
that the average number of atoms in a well nR^ ^ 1, but 
in the end we show that near the SI border this inequality 
follows from 

Let us first consider BCS phase corresponding to B > 
Bo- Here the criterion of superfluidity coincides with the 
condition that the Fermi level of weakly interacting Fermi 



nt 








Bn 



B 



FIG. 1: The phase diagram of the SI transition. Magnetic 
field B is plotted on the horizontal axis, while the fermion 
concentration n is plotted on the vertical one. S stands for 
superfluid and I for insulator. Bo is the Feshbach resonance 
point. Critical concentrations n/ and nb{B) for the classical 
random potential are straight lines given by Eqs. (jSj and (|10|) 
respectively. 
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gas reaches the mobihty edge in a given external poten- 
tial. Roughly speaking, this happens when the Fermi en- 
ergy of the gas Ep = (?i^/2m)(37r^n)^/'^ becomes larger 
than the amplitude Vq of the random potential. This con- 
dition leads to the critical concentration of SI transition 
n = Uf on the side of free fermions {B > Ba) 



CfR- 



Vo 



3/2 



(3) 



Thus, the segment of the SI border aX B > Bq is horizon- 
tal as shown in Fig.[T] It is known that in a classical long 
range potential the numerical coefficient C / can be found 
using the idea that in the classical potential ^ the mo- 
bility edge coincides with the classical percolation level 
Vp —1^. This is the level at which Fermi gas lakes formed 
in the random potential wells merge to create the infinite 
cluster or the Fermi sea. In a generic three-dimensional 
gaussian potential with the distribution function 



F{V) 



exp(-l^V2Fo2; 



this level corresponds to occupation of Oc 
space by lakes^. This gives 



(4) 



17% of the 



Vp = -0.96Vb. 



(5) 



Now we can find n f as the total concentration of fermions 
in wells deeper than Vp. Inequality ([2]) lets us use 
Thomas-Fermi (TF) approximation 



37r2 



{ 2m{Vp 



V) 



3/2 



F{V)dV (6) 



For a gaussian potential this leads to the coefficient in 
Eq. © 



Cf = 



37r5/2 



Vp/Vo 

[Vp/Vo-xf/^ exp{-x^/2)dx ~ 0.008. 

-oo 

(7) 

Let us switch to the less trivial BEC side of the diagram 
which corresponds to B < Bq. In this case, interaction 
of dimers plays the crucial role. Following Reffi we refer 
to the scattering length of the two dimers as add- Then 
the uniform gas of interacting dimers has the positive 
chemical potential 



4:TTh'^{n/2)add _ nh^nadd 
2m m 



(8) 



Here we took into account that the concentration of 
dimers is n/2, while the dimer mass is 2m. If fJ.{n) is 
larger than the amplitude of the random potential, Vq, 
the gas of dimers can screen the random potential redis- 
tributing a small fraction of its density from the hills of 
the random potential to the wells. On the other hand, if 
/i(n) <^ Vq the gas is fragmented in many disconnected 
lakes. Thus, the condition of delocalization of dimers 



and, therefore, the condition of superfluidity in this case 
is roughly speaking /i(n) = Vq. Substituting the nontriv- 
ial result of Ref.— 



add = 0.6a 



(9) 



into Eq. ([5]) and using Eq. ([T]) we get for the SI border 
concentration of fermions on the compact bosons side 
{B<Bo) 



miB) = CbR- 



,R Vq 
a fi^ /mR^ 



— ChR 



R 



Bn-B 



ao (h^/mR^) AB 
(10) 

Let us now estimate the numerical coefficient Cfc. To 
this end we again appeal to the percolation theory and 
deal with the percolation level 2Vp in the potential en- 
ergy of a dimer 2V(r). The local concentration n{r)/2 
of dimers adjusts to external potential according to the 
Gross-Pitaevskii equation (GPE) 



fiip{r) = 



4m 



2V{r) 



4:nh add 
2m 



(11) 

where is the condensate chemical potential, 2V{r) is 
the potential acting on a dimer, and the condensate wave 
function ^-){r) is normalized to total number of dimers, 
j dr\'4>{r)\'^ = N/2, where N is the total number of 
fermions. Thus, |7/>(r)p has the meaning of the local con- 
centration of dimers n{r) /2. Let us show that near the SI 
border one can use the TF approximation and drop the 
kinetic energy term of GPE. This can be done if the heal- 
ing length Ih = [{n/2)add]^^^^ of the condensate is much 
smaller than characteristic length of potential, Ih «C R. 
Using above estimate for the critical concentration ni, we 
get that for the classical disorder potential (Eq. ^) at 
the BEC side SI border 



Ihint) _ ( ImR 



R 



1^0 



1/2 



< 1. 



(12) 



Thus, one can proceed in the TF approximation, where 
at every point local concentration of the condensate n{r) 
satisfies equation 



•nfi nadd 



m 



+ 2V{r)=^l. 



(13) 



The chemical potential fi is determined by normalization 
of concentration of dimers n{r)/2 to the total number of 
dimers N/2 and grows with increasing N . li ji < 2Vp 
we get only disconnected Bose gas lakes. If /i > 2Vp the 
merging lakes form the Bose sea or the infinite cluster. 
Thus, similarly to the BCS side on the BES side the SI 
transition also happens when /i — 2Vp. For a gaussian 
potential with the help of Eq. ^ this gives for Ct in 

Eq. uni). 



Cb 



V2 



0.67r3/2 



[Vp/Vo - x] eiip{-x^/2)dx ~ 0.01. 

(14) 
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For more realistio^ distribution of the speckle potential 
we do not know the percolation threshold, but on the 
basis of approximate universality^ of the 9c we guess that 
Cb is the same as for gaussian potential within 20%. 

Thus, the SI border n{B) consists of the two straight 
lines, as shown in Fig. [TJ At _B < i?o it follows the 
line with the negative slope, Eq. (fTU]) and a,t B > Bq 
the border line is horizontal, Eq. ([3]). Eq. (jlOp is valid 
until nb{B) > Uf. At B ^ Bq - 5B, where 5B 

AB[Vb/(?iV™a§)]^^^ Eq. 1^ crosses over to Eq. ©. In 
the unitarity interval of the width SB around B = Bq at 
the SI border we arrive at n{B)a^ ^1. In other words 
in this interval kpa ^ Ih/ a ^ 1. Because the length of a 
dimer, ^ ~ add ^ ci, one can also say that only in the uni- 
tarity interval ~ 1 and dimers touch each other. This 
differs from the superconductor-insulator transition in a 
toy Coulomb model^, where the '--^ 1 at the long inter- 
mediate segment of the superconductor-insulator border. 

The fact that the critical concentration ni, exceeds nf 
is easy to understand. Indeed, at a given n dimers have 
much smaller chemical potential than weakly interacting 
fermions. Thus, dimers need a larger concentration n 
in order to get delocalized. In the similar way one can 
understand the growth of n{B) with the decreasing B 
at i? < Bq. Indeed, at a given n, the farther from the 
resonance, the more ideal the Bose gas of dimers is, the 
smaller is its chemical potential. Again, to compensate 
for this trend n{B) should grow with the decreasing B. 

Above we assumed that the number of particles in a 
well of the random potential is large, nR^ >> 1 and 
used the mean field approximation on the BEC side, ig- 
noring discreetness of particles. As we see from Fig.[l]the 
minimum value of the border concentration n{B) is nj. 
Therefore, inequality rifR^ ^ 1 guarantees that every- 
where on the border n{B)R^ ^ 1. Substituting Eq. ^ 
into UfR^ 1 we arrive at inequality Thus, this 
inequality is the single condition of validity of the above 
theory of SI border. 

It is clear from the above discussion that the insulating 
phase on the BEC side consists of disconnected lakes, 
populated by dimers. One can use the term Bose glass^™ 
for this phase, because it has no excitation gap and is 
compressible. 

Until now we assumed that the disorder is classical in 
the sense of inequality Let us now discuss what hap- 
pens for a quantum random potential, when the opposite 
strong inequality 



(15) 



takes place. In an experiment one can move from inequal- 
ity ^ to inequality (|15p by scaling down the intensity of 
the light beams creating speckles, while keeping the rest 
of the speckle set up (including R) fixed. How will then 
the phase diagram in {B,n) plane change? 

Let us start this discussion from the BCS side of the di- 
agram {B > Bq) and concentrate on the disorder induced 
density of states (DOS) at small energies. For simplicity. 



we assume that we are dealing with a gaussian poten- 
tial V{r) which two point correlation function decays as 

or faster at r i?. According to inequality 
the case of quantum disorder the wells of the size R do 
not have levels. In this case, the characteristic energy of 
the low energy tail of DOS is determined by wells of the 
size L ^ R, which are large enough to get a level^iiii^^. 
The typical depth of such wells V{L) is much smaller 
than Vo, namely V{L) = Vo{R/Lf/'^ < Vq. This hap- 
pens due to the cancellation of the majority of {L/R)^ 
contributions of wells the and hills of the size R. Using 
condition of the level existence V{L) = h^/mL^, we find 
the characteristic size of the well, which has a single level 



R 



iR^Vo 



(16) 



Substituting Lc to V{L) we arrive at the characteristic 
energy scale of the low energy tai fe^^i-^^ 



Vt = CtVo 



Vo 



h'/mR^ 



(17) 



The energy which separates localized and delocalized 
states (the mobility edge) is also of the order of Vt . For 
a quantum random potential (Eq. (fTSj) ) Vt <^Vo and the 
concentration of fermions which can be localized in the 
tails or, in other words, the critical concentration of SI 
transition, is very small, too 



R- 



Vo 



n^/mR'^ 



{Vo < h^/mR^). (18) 



Let us switch now to the BEC side of the phase diagram 
[B < Bo). In this case, the tails of DOS can accom- 
modate more dimers in the band of energies Vt because 
we can condense many bosons at one level. Only if the 
chemical potential of bosons given by Eq. ([5]) becomes 
larger than Vt the states become delocalized. Thus, rib 
can be estimated equating fi and Vt- Using Eq. ^ and 




FIG. 2: Schematic phase diagram of SI transition for a Fermi 
gas at B < Bo or for a weakly non-ideal Bose gas. The 
critical concentration of SI transition rib is plotted against 
the disorder amplitude Vo- S stands for superfluid and I for 
the insulator. The SI border shows crossover between regimes 
of quantum and classical random potentials at Vb = h^/niR^. 
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Eq. ini) we get for Vq < h^/mR'^ 

a\h^/mR'^J ao\h^/mR^J AB 

(19) 

This results can be also obtained from the conditioni^ 
that the SI transition happens at ^t/H ^ 1, where r 
is the relaxation time of a delocalized boson with the 
energy fi. The crossover between Eqs. ((T5| and (fT51) 
happens in the unitarity interval of the width SB — 
AB{ao/R)[Vo/{n^/mR'^)]'^ around Bq. In this interval 
na^ ~ fci?a ~ Ih/d ~ ^c/a ~ 1 similarly to the case 
of the classical potential. In the language of GPE the 
estimates we arrived above correspond to the solution, 
where all three terms in the right side of Eq. (fTT|) play 
comparable roles. In other words, expectations of the ki- 
netic energy term, of the random potential term and of 
the repulsion energy are of the same order of magnitude 
at SI border. Note, that at the same time the amplitude 
Vq of the bare potential is much larger than other terms. 
Only the quantum mechanical averaging makes the dis- 
order potential energy V{Lc) equal to other terms. For 
a Fermi gas, the idea of such averaging is known, for a 
long time^iiiii^. For a weakly non-ideal Bose gas idea 
of quantum mechanical "smoothing" of disorder poten- 
tial was explored only recently^i^^. However, SI phase 
diagram of an infinite, uniform in average gas could not 
be studied in Refs^ii^ because they dealt with an one- 
dimensional disorder potential. 

Let us discuss applicability of the mean field theory 
(GPE) for calculation of n?,. GPE is applicable if at 
n = ni, the characteristic length Lc or if nbL^ 3> 1. It 
is clear that rib^ rif. Multiplying this inequality by 
and using Eqs. and we arrive at necessary in- 
equality UbL^ 3> rifL^ ~ 1. Thus, the mean field theory 
is applicable for calculation of nb- Mean field approach 
fails and one arrives at single-particle regime^ only at 
n <C Ub- We see from Eqs. ^TE\\ and that in the case 
of a quantum random potential Vq <^ h /raR^ both crit- 
ical concentrations Uf and Ub decrease very rapidly with 



the decreasing Vq. As a result, while the whole phase 
diagram in this case looks like Fig. [TJ the concentrations 
Uf and nb are dramatically smaller than for a classical 
random potential. 

The summary of our results for the BEG phase is given 
in Fig. (21) where we plot the critical concentration of SI 
transition nb as a function of the amplitude of the ran- 
dom potential Vq (or the intensity of the speckle-building 
light), while scattering length a and the characteristic 
scale of disorder R are fixed. The fourth order parabola 
in the beginning of nf,(Vb) curve is given by Eq. ([TO]) . At 
Vo = Ti' jmR? this parabola crosses over to Eq. pH)) . Our 
results for the BEG phase of dimers shown on Fig. ([2]) 
are clearly applicable to a generic weakly non-ideal Bose 
gas with the scattering amplitude a > 0. 

In conclusion I have studied the zero temperature 
phase diagram of the superfluid-insulator phase transi- 
tion for a Fermi gas going through Feshbach resonance 
with the changing magnetic field and for a Bose gas. I 
dealt with uniform infinite gases, did not consider the role 
of the inverse parabolic profile n(r) in the trap, and did 
not study dynamics of the BEG phase expansion when 
the trap is eliminated. A likely scenario of this expansion 
(similar to that of Ref i^'^i^^ ) is as follows. If the maxi- 
mum concentration of the gas in the center of the trap 
TT-m < rib there is no expansion. On the other hand, if 
n„i > nb only bosons from the central domain r < tq 
where n(r) > nb leave reducing original n{r) to the flat 
n{r) = nb at r < Tq. This means that at n,„ ~ nb ^ nb 
the fraction of the gas mass released in an expansion is 
proportional to (nm — nbY^'^. 
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